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Abstract 

Stationary bifurcation control is studied under the assumption that the critical zero eigen¬ 
value is uncontrollable for the linearized system. The development facilitates explicit construc¬ 
tion of feedback control laws that render the bifurcation supercritical. Thus, the bifurcated 
equilibria in the controlled system are guaranteed stable. Both pitchfork bifurcation and tran- 
scritical bifurcation are addressed. The results obtained for pitchfork bifurcations apply to 
general nonlinear models smooth in the state and the control. For transcritical bifurcations, the 
results require the system to be affine in the control. 


1 Introduction 

In this paper, feedback control of stationary bifurcations is considered in the case that the critical 
mode is uncontrollable for the linearized system. The aim of the control design is described as follows 
Abed and Fu (1987). A nonlinear system is given for which a nominal equilibrium loses stability 
through a real eigenvalue crossing the imaginary axis at the origin. Under these circumstances, the 
system undergoes a stationary bifurcation. This can be a transcritical bifurcation or a pitchfork 
bifurcation. Figure 1(a) illustrates a transcritical bifurcation, and figure 1(b),(c) illustrate the 
two types of pitchfork bifurcation (subcritical and supercritical, respectively). The direction of a 
pitchfork bifurcation is its subcriticality or supercriticality. As discussed by Abed and Fu (1987), 
the supercritical pitchfork bifurcation (figure 1(c)) is preferable in practice since after the nominal 
solution has lost stability, new stable equilibria arise that provide a nearby operating condition. 
This observation motivates the search for feedback control laws that render a stationary bifurcation 
supercritical. This is the local stationary bifurcation control problem. 
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As shown in Abed and Fu (1987), nonlinear stabilizing controllers can be readily obtained if 
the zero eigenvalue (critical mode) is controllable for the linearized system at the bifurcation point. 
The situation was found to be significantly more difficult if the controllability condition fails, even 
though an analogous problem for Hopf bifurcation was addressed successfully (Abed and Fu 1986). 

Applications arise in which a stationary bifurcation occurs with the critical mode uncontrollable 
and for which the design of stabilizing controllers is of significant importance. One such example 
is control of rotating stall in axial flow compressors (Liaw and Abed 1996). 

There have been other investigations into design of bifurcation control laws for systems with an 
uncontrollable critical mode. For example, Fu and Abed (1993) consider the design of linear feed¬ 
back control laws for nonlinear systems affine in the control. Kang (1998) also studies stabilization 
for systems affine in the control under the condition that the noncritical modes are controllable. 

In the present paper, no assumption of controllability of the noncritical modes is made, and the 
general multi-input case is considered. Control of pitchfork bifurcations is considered for general 
nonlinear system models. For the case of transcritical bifurcations, the analysis is relegated to a 
narrow class of affine models. In both cases, the results permit explicit derivation of control laws. 

The remainder of the paper proceeds as follows. In Section 2, basic results on bifurcation analysis 
and control are recalled. In Section 3.1, these results are used for control design of general systems 
undergoing pitchfork bifurcation. In Section 3.2, control of a class of affine systems undergoing 
transcritical bifurcation is considered. Conclusions are collected in Section 4. 


2 Background 

This section reviews background material on stationary bifurcations and their control from Abed 
and Fu (1987). 

2.1 Bifurcation formulas for stationary bifurcation 

Determination of whether a stationary bifurcation is supercritical, subcritical or transcritical can be 
achieved using so-called bifurcation coefficients. These are coefficients in Taylor series expansions of 
quantities, especially eigenvalues, of bifurcated solutions in a small neighborhood of the bifurcation 
point. Formulas for these coefficients are referred to as bifurcation formulas. Next, we recall 
bifurcation formulas that will be needed in the control design of this paper. References (Howard 
1979, Iooss and Joseph 1980, Abed and Fu 1987) can be consulted for further details. 
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Consider a one-parameter family of nonlinear autonomous systems 


x = f(x,n) (1) 

where x 6 R n is the state vector and fi is a real-valued parameter. Let /(x,/t) be sufficiently 
smooth in x and /j and let x ® be the nominal equilibrium point of the system as a function of the 
parameter fi. 

Suppose that the next hypothesis holds. 

(S) The Jacobian matrix of system (1) at the equilibrium x° has a simple eigenvalue Ai (//) such 
that Ai(0) = 0 and A^O) ^ 0, and the remaining eigenvalues lie in the open left half of the 
complex plane for /t = 0. 

The stationary bifurcation theorem (Chow and Hale 1982, Guckenheimer and Holmes 1983) 
asserts that hypothesis (S) implies a stationary bifurcation from x® at p = 0 for (1). That is, a new 
equilibrium point bifurcates from Xq at p = 0. Near the point (xq,0) of the (n + l)-dimensional 
(x, p)-space, there is a parameter e and a locally unique curve of critical points (x(e), //(e)), distinct 
from xj) and passing through (xq,0), such that for all sufficiently small |e|, x(e) is an equilibrium 
point of (1) when p = //(e). 

The parameter e may be chosen such that x(e),//(e) are smooth. The series expansions of 
x(e),/t(e) can be written as 


M( e ) — Mi e + /^2e 2 + • • • (2) 

x(e) = x° + xie + x 2 e 2 + ... (3) 

If fii ^ 0, the system undergoes a transcritical bifurcation from x° at // = 0. That is, there is a 
second equilibrium point besides x° for both positive and negative values of p with |/x| small (see 
figure 1 (a)). The stability of the bifurcated equilibria in the case of transcritical bifurcation is as 
depicted in the figure. If = 0 and )i 2 ^ 0, the system undergoes a pitchfork bifurcation for |/t| 
sufficiently small. That is, there are two new equilibrium points existing simultaneously, either for 
positive or for negative values of /t with //. small (see figure 1 (b),(c)). The bifurcated equilibrium 
points have an eigenvalue /3(e) determining their stability which vanishes at e = 0. The series 
expansion (3(e) is given by 

(3(e) = (3\e + foe 1 + ... (4) 

with 

fo = -/hA'(0) (5) 
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and, in case d\ = 0 , /3 2 is given by 

fh = -2/x 2 A'(0) (6) 

The stability coefficients (3\ and /3 2 can be determined solely from eigenvector computations 
and the coefficients of the series expansion of the vector held. System (1) can be written in the 
series form 

x = L^x + Q^(x,x) + C^(x,x,x) + ... 

= Lqx + pL\x + ^ 2 L 2 x + ... 

+Qo(x,x) + nQ±(x,x) + ... 

+C 0 (x,x,x) + ... (7) 

where x = x — Xq,L m ,Li,L 2 are n x n matrices, Qfj,(x, x), Qo(x, x), Qi(x, x) are vector-valued 
quadratic forms generated by symmetric bilinear forms, and C fJj (x,x,x),Co(x,x,x) are vector¬ 
valued cubic forms generated by symmetric trilinear forms. 

By hypothesis (S), the Jacobian matrix Lq has only one simple zero eigenvalue with the re¬ 
maining eigenvalues stable. Denote by l and r the left (row) and right (column) eigenvectors of the 
matrix Lq associated with the simple zero eigenvalue, respectively, where first component of r is 
set to be 1 and the left eigenvector l is chosen such that Ir = 1. It is well known that 

A'(0) = lL\r ( 8 ) 

A stability result for the bifurcated equilibria of system (7) is given in the following lemma. 

Lemma 1 If [3 1 = 0 and /J 2 7 ^ 0, then Eq. (1) undergoes a pitchfork bifurcation at n = 0. ///?i 7 ^ 0, 
then Eq. (1) undergoes a transcritical bifurcation. In the former case, the pitchfork bifurcation is 
supercritical if @2 < 0, but is subcritical if @2 > 0. Here, 


fh 

= lQo(r,r) 

(9) 


= 2l{2Q 0 (r,x 2 ) + Co(r,r,r)} 

(10) 


where 


X2 


-(R t K)- 1 R t 


Qo{r,r) 

0 
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( 11 ) 


Lo 

l 





(a) Transcritical bifurcation (b) Subcritical pitchfork bifurcation 


-:stable 

-runstable 

(c) Supercritical pitchfork bifurcation 

Figure 1: Stationary bifurcation diagrams. The variable in the vertical direction represents equi¬ 
librium values for a system in normal form. 



2.2 Bifurcation control in the case of a controllable critical mode 

Next, we recall the results of Abed and Fu (1987) which give sufficient conditions for local stabi- 
lizability of an equilibrium point at criticality and for local stabilizability of bifurcated equilibria 
under assumption (S). These conditions involve assumptions on the controllability of the critical 
mode of the linearized system. 
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Consider a system 

x = f(x,n,u) (12) 

where x € R n is the state vector, u € R m is input vector, and fi £ R is the bifurcation parameter. 
Expanding the vector field of (12) in x, // and u, we have 

* = f{x,n,u) 

m m 

= L 0 x + jiLix + y^,UjL\x + b l Uj + Qo(x,x) 
i=l i=l 


m 


+^ 2 L 2 x + nQi(x, x) + ^2 UiQi(x, x) 

1=1 


+Cq{x, X, x) + ... 

(13) 

Note that the linearized system is 


m 

x = Lqx + ^2 b l Ui 

(14) 


i=l 


At the critical parameter /j := /j c , Lq has one zero eigenvalue. A feedback control consisting of 
quadratic and cubic terms is considered: 

m 

u{x) = ^2x T Q Ui x + C Ui (x,x,x) (15) 

i =1 

Following (Abed and Fu 1987) use an asterisk to denote quantities for the closed-loop system. Then 
the bifurcation coefficients taking into account the state feedback (15) are (details in (Abed and 
Fu 1987)): 


Pi = l{Qo(r,r)+J2( rT Qui r )b 1 } (16) 

i= 1 


and, if 8f = 0, 


P* 2 = 2l{2Q 0 (r,x* 2 ) + 2j2(r T Qu i x* 2 )b i 

i= 1 


m m 

+ ^ ^2( rT Qui r L\r ) + C 0 (r, r,r) +J2 C n,i r , r, r)tf} 

i =1 i=l 


Here, x 2 is given by 


x 2 


-{R t R)- 1 R t 


Qo(r,r) + E r r=iir T Q Ui r)R 

0 


(17) 
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( 18 ) 


m 

= x 2 -Y,(r T Q u . i r)(R T R)- 1 R T 

i=1 

where x 2 is the value for u = 0. 

By employing these formulas, the following two results were obtained in (Abed and Fu 1987). 

Theorem 1 Let hypothesis (S) hold and assume lb 1 ^ 0 for some i € 1 ,m, that is, the critical 
zero eigenvalue is controllable for the linearized version of Eq. (12) near the origin. Then there is 
a smooth feedback control u = u(x) with u( 0 ) = 0 , containing only quadratic and cubic terms in 
x, which solves the local stationary bifurcation problem for (12). Moreover, the quadratic terms in 
u(x) can be used to ensure that j3\ = 0 for the controlled system, and the cubic terms can then be 
used to ensure that f3 2 < 0 . 

Theorem 2 Let hypothesis (S) hold and assume lb 1 = 0 for all i = 1 ,m, that is, the critical 
zero eigenvalue is uncontrollable for the linearized version of Eq. (12) near the origin. Then if 
(3\ ^ 0 for with u(x) = 0, the local stationary problem for (12) is not solvable by a smooth feedback 
control with vanishing linear part. 


V 

0 


3 Stationary Bifurcation Control with Uncontrollable 
Linearization 

In this section, we investigate the design of feedback control laws for stabilizing stationary bifur¬ 
cation for systems with an uncontrollable critical mode. We first consider the case of pitchfork 
bifurcation and then transcritical bifurcation. 

3.1 Pitchfork bifurcation control 

Suppose that the system (12) undergoes a pitchfork bifurcation and that the critical (zero) eigen¬ 
value is uncontrollable for the linearized system at bifurcation. This implies that the bifurcation 
coefficient (3 1 = 0 for the open-loop system. 
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To facilitate the search for a stabilizing control law, perform a linear coordinate transformation 


z 


Tx 




(19) 


where r, £ 2 ,... are orthogonal column vectors in R n . Thus det T ^ 0. In z coordinates, the 
system becomes 


z = 


m m 

= L 0 z + fxLiz + ^2 UiL\z + ^2 b l Ui + Qo(z, z) 

i=l i=1 


m 

+H 2 L 2 z + nQi(z,z) + ^2uiQi(z,z) 

i =1 


+Co {z, z, z) + ... 


( 20 ) 


For simplicity, we do not change the notation for the quantities Lo, L\,Qq, ..even though they 
are affected by the transformation. Because of the orthogonality assumption above, L$ takes the 
form 


To — 


0 0 
0 Ai 


( 21 ) 


where A\ is an (n — 1) x (n — 1) invertible matrix. Take the right and left eigenvectors corresponding 
to the zero eigenvalue of Lo (21) to be r = [1,0 ,..., 0] T and l = [1,0,..., 0 ], respectively. 

Since Pi = 0 in the absence of control, we limit our search to feedback laws with a vanishing 
linear part. This ensures that /3J, the value of Pi after feedback, is also 0. (As in the preceding 
section, an asterisk indicates a quantity for the closed-loop system.) Thus, controls are sought in 
the form 

Ui{z) = z T Q Ui z + C Ui (z,z,z) for i = 1,..., m (22) 


The coefficient P 2 , which determines whether the pitchfork bifurcation is supercritical or subcritical 
for the closed loop system, is given by (17). Using l = [1, 0,... , 0] and the structure of Lq, we find 


that 


{R T R)~ l R T 


0 0 1 
0 Af 1 0 


(23) 


In addition, since only the first element of r is nonzero, 


r T Q Ui r = c l u 


(24) 
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where c l u denotes the coefficient of z\ in the formula (22) for the input «j. Thus, we can conclude 
that only the z\ terms in an input of the form ( 22 ) affect the direction of the pitchfork bifurcation. 
Eq. (18) for x\ now gives 

m 

x *2 = x 2 ( 25 ) 
i=1 

where the vectors i = 1,..., m, are defined as 


1 

0 i 

1 _ 


0 

KJ 




(26) 


and where the b l are (n — l)-dimensional vectors containing all but the first ( 0 ) component of the 
vectors 6 *; i.e., 


¥ = 


0 

¥ 


(27) 


Note that lb 1 = 0 and l = [1,0,..., 0]. 

Again, using the fact lb 1 = 0 for i = 1,..., m, Eq. (17) simplifies to 


ft = 2l{2Q 0 (r,x* 2 ) + Y,ftL\r + C 0 (r,r,r)} (28) 

i=1 


Note that the eigenvectors l and r are unaffected by the feedback since the control law does not 
contain linear terms; thus the asterisk notation is not needed for these vectors. This will change 
below in our discussion of control of transcritical bifurcations. Note also that the cubic terms in the 
input do not have any effect on ft. Moreover, since only the first elements of r and l are nonzero 
and the z\ term in z\ vanishes (since ft = lQo(r,r) = 0), Eq. (28) gives 


n m n m 

& = 4{£ Cux \2 - £4£ c i iHi + 2 (£4^ i 1 ) + 2c 3 
i =2 j =1 i —2 j =1 

Tl 7TL TL 

= 2{2 £ cux l 2 + 2c-j } + 2 £ ft {£ c\jCi + ft 1 } ( 29 ) 

i =2 i= 1 j =2 


where cu and C 3 are the coefficient of z\Zi and the coefficient of z\ in z\, respectively. Also, L ) 1 
denotes the (1,1) element of the matrix L\. 

The following assumption is now introduced. 


9 



(Al) At least one among the coefficients p*, i € 1,..., m, does not vanish, where the pi are defined 
as follows: 

n 

Pi = -4^ Clj Cj + 2L { 1 (30) 

3 =2 

Theorem 3 Let the system (12) with input u = 0 undergo a pitchfork bifurcation from the origin at 
p = 0 . Also, assume lb 1 = 0 for alii = 1 ,..., m; that is, the critical zero eigenvalue is uncontrollable 
for the linearized version of (12). Moreover, suppose that (Al) holds. Then, there exists a smooth 
feedback control containing only quadratic terms in z\, which solves the local stationary bifurcation 
control problem for Eq. (12). Moreover, the nonvanishing of at least one of the pt is necessary and 
sufficient for stabilization using feedback of the form (22). 

3.2 Transcritical bifurcation control 

Next, we consider the case in which the system (12) undergoes a transcritical bifurcation, that is, 
(S) holds and /?i 7 ^ 0. As in the foregoing, the critical zero eigenvalue is assumed uncontrollable for 
the linearized version of Eq. (12). We proceed in two steps. First, we use linear state feedback to 
transform the transcritical bifurcation into a pitchfork bifurcation. Then, we add quadratic (plus 
higher order if desired) feedback to stabilize the achieved pitchfork bifurcation. 

For simplicity, we again consider the system after the linear transformation (19) has been 
performed. From Theorem 2, a feedback consisting of only quadratic and cubic terms cannot 
transform the transcritical bifurcation into a pitchfork bifurcation. It is for this reason that we seek 
a feedback control that also includes a linear term, taking the form 

Ui(z) = KiZ + z T Q Ui z + C Ui (z, z, z) for i = 1,... ,m (31) 

Unfortunately, a linear term in the state feedback significantly complicates analysis of the bifur¬ 
cation. A linear feedback can affect the eigenvalues and eigenvectors of the linearization, making 
determination of its effect on the bifurcation formulas (9), (10) very difficult. 

To mitigate this problem, we focus on the following special class of nonlinear systems: 

m 

z = f(z,p) + ^b'ui (32) 

i=1 

Note that this model is affine in the inputs, and that the inputs enter the dynamics through 
multiplication by constant vectors. 
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For the system (32) to undergo a transcritical bifurcation for open-loop system, bifurcation 
coefficient (5\ must be nonzero. The bifurcation coefficient fi\ for the open-loop system (32) can be 
calculated by using Eq. (9). Recall that only the first element of the left and right eigenvectors 
corresponding to critical zero eigenvalue are nonzero. The coefficient (3i is given by 


ft 


lQo{r,r) 


1 0 ... 0 


T— 

r r 

T — 1 

r 

T— 

r ^ n r 



(33) 


where the Hj are symmetric matrices that are easy to compute given Q o (the bilinear form Qo(x,y)). 

Since the first element of each of the b l vanishes (27), the controls m do not change the right 
side of the equation for z\ in (32). Thus, Hi is not affected by feedback. It should also be noted 
that linear state feedback does not change l (the left eigenvector associated with the critical zero 
eigenvalue) due to the structure of the Jacobian matrix (21) and the uncontrollability of the critical 
mode. Thus, the only way to transform the transcritical bifurcation into a pitchfork bifurcation 
(i.e., to make (3\ = 0) is to change the right eigenvector r through linear state feedback. 

Here, for simplicity we limit the linear state feedback to be a function of z\ alone. Thus, we 
consider state feedbacks of the form 


Ui(z) = kiZ\ for i = 1,..., m (34) 

where ki G R for i = 1 ,,m. (Note that the i-th scalar k{ is the first component of the row vector 
Ki of Eq. (31), for i = 1,... , to.) Once appropriate gains ki are found, nonlinear terms will be 
re-inserted in the feedback. 

In Theorem 4 below, we will require at least one of the following two assumptions to hold. 


(Bl) There exists a v 


Vl 

V 


with v\ 7 ^ 0 vi € R, that belongs to the null space of Q i and v 


belongs to the linear space spanned by the for i = 1 ,..., m (26). 


(B2) Q i has a positive eigenvalue (A p ) and the negative eigenvalue (denote A n ) with corresponding 
eigenvector p and n such that following three conditions hold. 


(B2-1) At least one of the eigenvectors {n and p) has a nonzero first element. 


(B2-2) Writing n = 


ni 

and p = 

Pi 

n 


P 


where n\ and p\ are scalars, n are p belong to the 


linear space spanned by the Q l for i = 1 ,..., m (26). 
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(B2-3) The following pair of equations is valid: 

ni+pi = 1 
n n 

X pJ2p 2 i + = 0 ( 35 ) 

i =1 i=l 

where p* and rii denote the i-th element of p and n, respectively. 

The condition (B2-3) fixes a scaling on the eigenvectors n and p. Also, note that if only one of n\ 
and pi is nonzero, then Eq. (35) can always be achieved. 


Theorem 4 Let the system (32) with input u = 0 undergo a transcritical bifurcation from the 
origin at p = 0. Also, assume lb 1 = 0 for all i = 1 ,... ,m; that is, the critical zero eigenvalue is 
uncontrollable for the linearized version of (32). Moreover, suppose that either assumption (Bl) 
or (B2) holds. Then, the transcritical bifurcation of system (32) is transformable to a pitchfork 
bifurcation by means of linear state feedback of the form (34). 

Proof: The Jacobian matrix after applying a linear state feedback of the form (34) is given by 


L 


* _ 

0 — 


o o 

E™ i hCi Ai 


(36) 


The right eigenvector r corresponding to the zero eigenvalue becomes 


where the first element 


of r is 1 since Ir = 1 and l = [1, 0 ,..., 0], and where the subvector r solves the following equation: 


53 W + M? = 0 (37) 

i =1 

Since A\ is an (n — 1) x (n — 1) full rank matrix, we have 

m 

r = 53 kiAr 1 # 

i =1 
m 

= Y, k iC ( 38 ) 

i =l 

It is easy to see that if either assumption (Bl) or (B2) is satisfied, then we can set the right eigen¬ 
vector corresponding to the critical eigenvalue to a value which results in (3\ = 0. ■ 
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Note that including other states besides z\ in the linear state feedback will change the A\ 
matrix. Thus, it is possible to design a linear state feedback that will maintain stability of matrix 
A\ and result in a desired set of vectors Cf = A^ 1 b l which satisfies the conditions of Theorem 4. 
The following theorem follows directly from these observations. 

Theorem 5 Let the system (32) with input u = 0 undergo a transcritical bifurcation from the origin 

at p = 0. Assume lb 1 = 0 for all i = 1 ,..., m; that is, the critical zero eigenvalue is uncontrollable 

for the linearized version of (32). Moreover, suppose that Si is either positive definite or negative 

definite. Then, bifurcation control of the system (32) cannot be achieved by means of state feedback. 

Proof: For a pitchfork bifurcation to occur in the closed-loop system, (5\[ must vanish. For the 
system (32), is given by 

P{ = r* T Sir* (39) 

where r* is the right eigenvector of Lq (the closed loop Jacobian matrix) associated with the eigen¬ 
value 0 . Since Si is either positive definite or negative definite, there is no nonzero r* G R n for 
which PI vanishes. ■ 

Now, we want to add a quadratic term and/or cubic term to the linear feedback such that the 
pitchfork bifurcation that system (32) will undergo is guaranteed to be supercritical. To determine 
direction of pitchfork bifurcation, we calculate p 2 . Recall that the left eigenvector corresponding 
to the critical zero eigenvalue is unaffected by linear state feedback. Also, recall that only the first 
element of the left eigenvector is nonzero and it takes the value 1. Using these facts and Eq. (17), 
we find 

m 

P% = 4r* T E 1 x* 2 + 2'£r* T Q Ui r*L\ i r* + 2C 01 (r*,r*,r*) (40) 

i= 1 

where C/u represents the cubic terms in z\ and L l p denotes first row of matrix L\. Using Eq. (18), 
(40) becomes 

m m 

P* 2 = M* T Z 1 x 2 -4r* T Z 1 '£r* T Q Ui r*5 i + 2'£r T Q Ui rL\ i r* + 2C 01 (r*,r*,r*) 

i= 1 i= 1 

m 

= M* T Z 1 x 2 + 2C 01 (r*, r *,r*) + J2 r * T Qu i r *{-M* T Z 1 6 i + 2L¥r*} (41) 

i =1 
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where 8 t = {R T R) l R T 


V 

0 


The next theorem follows from the preceding discussion. The 


theorem make use of the following assumption. 

(B3) At least one among the coefficients i G 1,... ,m does not vanish, where the i>i are defined 
as follows: 

Vi := -4r T £i^ + 2(42) 


Theorem 6 Let the system (32) with u = 0 undergo a transcritical bifurcation from the origin at 
H = 0. Assume lb 1 = 0 for all i € 1 ,... ,m, that is, the critical zero eigenvalue is uncontrollable for 
the linearized version of (32). Moreover, suppose that (B3) holds and also that at least one of (Bl) 
or (B2) holds. Then, there exists a smooth feedback control in the form of (31), which solves the 
local stationary bifurcation control problem for Eq. (32). 

Proof: Let i be an index for which v-- ^ 0. Then, letting the control u-■ of (31) be u^(z) = 
K^z + z t Q u ,z = kjZi + z t Q u ,z we can set /?| (41) to any desired value. One such example, which 
is independent of r*, is z t Q Uo z = cz\ since the first element of r* is always 1. Here, c is a real 
constant chosen to ensure that /3| is negative. ■ 


4 Conclusion 

In this paper, we studied stationary bifurcation control for systems with uncontrollable lineariza¬ 
tion. For systems undergoing pitchfork bifurcation, purely nonlinear state feedback was used to 
achieve a supercritical bifurcation in the closed-loop system. This was done under general conditions 
for a large class of real-analytic models. Analogous results for control of transcritical bifurcation 
were obtained for a smaller class of systems using a two-step approach. In the first step, linear 
feedback was used to transform the bifurcation into a pitchfork bifurcation. In the second step, 
quadratic feedback was employed to ensure supercriticality of the achieved pitchfork bifurcation. 
In both cases, the designed feedback laws would represent the first phase of an actual design, which 
could be followed with inclusion of other feedback terms to optimize system performance. 
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